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. We consider the origin of new phases in super symmetric grand unification 

model, and show how significant new contributions arise from the gluino me- 
diated diagram. We then present a more general model independent analysis 
of various modes of B-decays suggested previously for measurement of the 
CKM phases and point out what they really measure. It is in principle pos- 
sible to separate out all the phases. 
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We consider the origin of new CP violating phases from physics beyond the standard 
model (SM) and their effect on various measurements of CKM phases a, (3 and 7 proposed 
hitherto |0]-f§|. Among new sources of CP violation are multi-Higgs models ||, the left- 
right model U and supersymmetry. In this note we focus on supersymmetry, which is 
very attractive from a grand unification viewpoint and provides many new sources of CP 
violation. One obvious source is the complex soft terms. Even when these are taken to be 
real, unification of right handed fields, like the left handed ones, can lead to a new source 
of CP violation. For example, a group like SO (10) f^H or models with intermediate gauge 
groups 0,0 like SU(2) L x SU(2) R x SI/(4) C , SU(2) L x SU(2) R x SU(3) C x U(1) B -l have 
these extra phases. Supersymmetric contributions with new phases can be as large as the 
SM in the B — B mixing and loop processes that lead to b — > sqq. 

In this paper we make the first complete calculation of the gluino contribution to Amg in 
a SUSY grand unified S0(10) theory. This calculation can easily be extended to the models 
with the intermediate gauge symmetry breaking scale considered in references |10| , |TT| . This 
calculation has been done previously by assuming same masses for the SUSY particles only 
for the low tan/5 scenario ||. We consider two scenarios: (i) where the Yukawa couplings 
are unified (i.e. large tan/3 scenario) and (ii) a low tan/? scenario. We show how the new 
contributions are large and can affect the interpretation of measurement of CKM phases. 
We then discuss the specific B-decay modes needed to extract the CKM phases even in the 
presence of new physics. This discussion actually uses model independent analysis that is 
valid in almost any kind of departure from the SM. 

Since the soft SUSY breaking terms are gravity induced, we shall assume them to be 
universal at the scale 2.4 • 10 18 GeV(M x ) which is the reduced Plank scale. For simplicity we 
also assume the soft terms to be real. It has been shown that a grand unified model based 
on SO(10), which we will use in this paper, gives rise to flavor violating processes in both 
quark and lepton sectors. Consequently lepton flavor violating processes like /1 — > put 
bounds on the parameter space along with b — > 37 [IB|. For models with the intermediate 



gauge symmetry breaking scales, the soft terms can be universal even at the GUT scale and 



still give rise to these effects. The superpotential for the Yukawa sector at the weak scale 
for the SO(10) grand unification or for the grand unifying model with an intermediate scale 
can be written as |T|-|9|,|rT| : 



W = Q\ U XJ C H 2 + QV*A d S 2 V t £> c ffi + £ c V^A l S 2 v2jL#i , (1) 

where V is the CKM matrix, is the CKM matrix at the GUT scale (for intermediate 
gauge symmetry breaking models G is replaced by I to denote the intermediate scale) and S 
is the diagonal phase matrix with two independent phases. The phases in the right handed 
mixing matrix for the down type quarks and the down type squarks can give rise to new 
phases in Am# and Am^ through the gluino contribution. 

The existing calculation [14j] for Am# using GUT model usually assumes that the soft 



terms are universal at the GUT scale(~ 10 16 GeV). Under that assumption it is found that 
charged Higgs has the dominant contribution. But with the universal boundary condition 
taken at the Planck or string scale there can be a large contribution from the gluino mediated 
diagram due to the fact that the fields that belong to the third generation have different 
masses compared to the other generation at the GUT scale due to the effect of the large top 
Yukawa coupling which gives rise to the non-trivial CKM like mixing matrix in the right 
handed sector. We first consider large tan/3 solution. In order to have a realistic fermion 
spectrum and the mixing parameters in the large tan /3 case, we use a maximally predictive 
texture developed in the reference |I5|| . We will look at a scenario where X^Mq) = 1 and 
tan/3 = 57.15, which gives m 4 =182 GeV and mft=4.43 GeV. For the small tan/5 scenario we 
have used Xt(Ma) = 1-25 and tan/? = 2. Above the GUT scale we use one loop RGEs for 
the soft terms and the Yukawa couplings [§] . Below the GUT scale we will use the one loop 
RGEs in matrix form in the 3x3 generation space for the Yukawa couplings and soft SUSY 



breaking parameters as found in Ref. p4| , |T6|| rather than just running the eigenvalues of these 
matrices as is often done. Although doing this does not provide any new information when 
tan p is small, when tan (5 is large it allows one to know the relative rotation of squarks to 
quarks and sleptons to leptons. In the large tan (3 scenario to make sure that the electroweak 



symmetry is broken radiatively, we need a non zero value of the D-term usually referred as 
m 2 D , which gets introduced at the GUT scale due to the reduction of rank in the SO(10) f|T7|] . 
This D-term is also bounded from above and below by requiring the pseudo scalar mass to 
be positive along with the squark and slepton masses ||13|| . 

We calculate Am# using gluino contribution and compare it with the SM result. We 
have done the calculation in SO(10), though this calculation can easily be generalized to the 
models with the intermediate scale and other grand unifying models. We use the expression 
for ArriK given in the reference [JTSj] (modified for the purpose of B° — B° mixing), because 
these expressions use the squark mass eigenstate basis derived from the full 6x6 mass 
matrices, which automatically incorporates mixing between "right-handed" down squarks 
and right-handed down quarks as is inevitable with either large tan/3 or SO (10) grand 
unification. We plot the ratio ArriB gluino / AmB SM as a function of \x for different values of 
the gaugino masses (miw) in Figure 1 for the large tan/3 case, where mo is 1 TeV for the 
entire plot. The gluino mass is related to the gaugino mass by the relation rrig = We 
take a s (M z ) = 0.121 and ac = 1/23.9 and the scale for grand unification to be Mq = 2 ■ 10 16 
GeV. Also in this scenario we have three variables: m (the universal scalar mass), m 1( / 2 (the 
universal gaugino mass) and the m 2 D (throughout our analysis we will assume the trilinear 
soft SUSY breaking scalar coupling A = at the Planck scale). The upper and the lower 
end of each curve correspond to the upper and the lower limit of the D-term respectively. 
As mentioned in the reference ||13|| , the parameter space with m less than 1 TeV as well as 



fi > is restricted by the flavor changing neutral currents. In Figure 2 (small tan (3 case) we 
plot r(= Am b luino / Arris SM ) as a function of the gaugino mass (m^) for different values of 
the scalar masses mo, where tan/? is assumed to be 2 and \x < 0. In the plot we have used 
the absolute value of Am^. We restrict ourselves to the parameter space allowed by the 
other flavor changing decays. We also make sure that /i is less than 800 GeV to avoid fine 
tuning. In both figures the SUSY contribution can be comparable to the SM. As a matter 
of fact in this parameter space the gluino contribution to the b — > S7 is also large P~2"| , |T3"| . 



From the graph one can see that for the scalar mass (or the right handed slepton mass) 



mo=l TeV and for the gaugino mass mi > 140 GeV (or the gluino mass > 405 GeV), the 
SUSY contribution is small (less than 20 %) compared to the SM in the large tan (3 scenario. 
In the low tan/5 scenario, for the gaugino mass mi > 200 GeV (or the gluino mass > 578 
GeV) and for the scalar mass (or the right handed slepton mass) m > 200 GeV, the SUSY 
contribution becomes small (less than 20 %) compared to the SM. For a complete SUSY 
calculation, there could be contributions from charged Higgs, chargino and neutralino. The 
charged Higgs contribution does not change significantly with the new boundary condition 
and has been found to be comparable to or even greater than the SM contribution when 



the soft SUSY breaking terms are taken at the GUT scale jlj]. Also this contribution does 
not involve any right handed down type quark-squark mixing, so that it has the same phase 
structure as the SM does. Consequently, the CKM measurement is not affected from the 
charged Higgs contribution as we will discuss later. Chargino and neutralino contributions 
are usually small JT4]JT9|] and have no effect on the CKM measurements. 

The soft terms (e.g A and or fi) can also be complex. In that case one can get phases in 
ArriB even without grand unification. The complex terms in the mass matrix for the squarks 
and sleptons are then responsible for the new phases which are somewhat restricted by the 
edm of electron or neutron |^D[ , however large phases can appear when the scalar masses are 



in the TeV range [pijl . There could also be an induced phase in A due to the phase in the 
Yukawa sector through renormalization, even when A is real at the GUT scale. The phase 
induced is really small and gives rise to the edm of electron well within the experimental 
limit for squarks and gluino masses O(100 GeV) ||22|| . It is possible to get comparable Am^ 
and Arrix from supersymmetric contribution with new phases |23[ in a model based on the 



MSSM (without grand unification) with right handed mixing matrix in the up sector. 
The contribution to the Am^ s can be parameterized as : 

Am Bd = Asm + B SU sy + Cm/sre*, (2) 
Am Bd =A Bd V t iv t 2 b e^. 

To make the analysis a most general one we have included Bsusy which has the same phase 



structure as Asm- in our example, for Am^, the box diagram with the LRLR structure 
(helicities of the fermions in the external legs) has the mixing structure \V td \ 2 e 1 ^, and RRRR 
type of box diagram has the mixing structure V^e 2 ^ in the diagonal quark mass basis with 
just b squark in the loop, where arises from the matrix S. Note that even if <ft is 0, both 
the RRRR type and LRLR type still have different phase structures compared to the SM. 
As a matter of fact any contribution from beyond the SM including multi-Higgs models 
and left-right models can be written as above. AB d e l ^ B d originates from the combination of 
the SM contribution and the new contribution. Similarly we have for B s — B s and K — K 
mixing: 

Am Bs = AbXsV^* , Am K = A K v£v&**. (3) 
Expressions for q/p for each of these mesons are now: 

\p) Bd \v»vz) ' \p) Bs {v tb v t *J ' \p) K \v cd v*J 1 ; 

In general <pB d i 4>B a and 4>k are unrelated. These phases are so defined that they are in ad- 
dition to the phases present in the SM, and can be treated as separate observables. Charged 
Higgs mediated box diagram has C = 0, and CKM measurements are unaltered. However, 
in our calculation C is non-zero (the LRLR type and the RRRR type), and will affect CKM 
phase measurements. 

The CKM phases are defined as: a = Axg{-V* b V td /V* h V ud ) , (3 = Arg(-V; b V cd /V t * b V td ) , 
7 = Arg (—V* b V ud /V* b V cd ). Based on the SM, many methods have been suggested for measur- 
ing these CKM phases using B decays |]-§J- The cleanest method involves time- dependent 
measurements of rate asymmetries in neutral B decays to CP eigenstates 0] , where one mea- 
sures the time-dependent rate asymmetry a,f cp (t) which is a function of A = (q/p) i (A/A^j, 
where i denotes the corresponding mixing, i.e., i = B d , B s , or K, and A = A(B° — ► /) and 
A = A(B° -> /) with CP eigenstate /. 

We shall analyze the different CP eigenstates that have been suggested, and consider 
carefully what phases the measurements now yield. Our assumption for decay amplitudes 
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is that, while the tree amplitudes have the SM phases, any loop process could have an 
additional unknown phase arising from beyond the Standard Model. Thus for penguin 
amplitudes we have 

— — ( —] p i<Ppeng (r\ 

^ \A) SM e ' (5) 
where 4> pen g is a phase in addition to SM phase. The results of our analysis are presented 
in a convenient tabular form (Table I) modeled after a similar table in the SM given in 



24j. Some of the modes have also been discussed in the reference || where only the SUSY 
grandunification contributions are retained. It is important to realize that with our defi- 
nitions of additional phases as defined in Eqs. (4) and (6), these phases are measurable. 
Further, the analysis is essentially model independent, as these new phases can arise in any 
model beyond the SM. 

In row (1) we consider B d — > ipKs- This mode which is tree dominated has ImX given 

by: 

V t b v td\ e i4, Bd { Vc S Vcd \ e i<t, K f V cbV* Q 



ImX = Im 



Im 



v tb v td j \v*v cd j \v; b v CS/ 

= -sin(2/3 + B(i + <^). (6) 

Note that the mode b — > ccs has negligible penguin contribution. In the SM this measure- 
ment yields sin (2/5). Similarly B s — > ifxp, iprj would yield <pB s while in the SM there is no 
asymmetry. In row (2) and (4) we have pure penguin processes b sss and b — > dds, 
respectively. These could have an additional weak phase (p pe ng or 4>' peng corresponding to 
each process. In row (2) the weak phases in B s and Bd are the same 4> peng because they arise 
from the same quark subprocess. The processes in row (3) are generally not suitable as both 
tree and penguin amplitudes make comparable contributions to the final states. In row (5) 
tree amplitude dominates and although the modes are Cabibbo suppressed, they are useful. 
In row (6) it is assumed that in the SM top contribution dominates in the loop. The contri- 
butions from charm and up quarks are expected to be about 10% over most of the allowed 



range p5 |. In row (7) tree contribution dominates and the small penguin admixture can be 



removed using isospin analysis [§]. Row (8) has processes dominated by tree diagrams and 
even though the mode D°K* is not a CP eigenstate, an analysis of this mode can be used 
to determine 7 ||. The charged B decay mode D°K + can be used alternatively, based on 
the same type of analysis 

It is clear from the table I that from Bd decays we can extract the combination (3 + cj)B d /2 
and 4>k, 4>peng and 7. From B s decays it is possible to measure 4>b 3 , (fix, (fipeng, 4>'peng an d 7 
and the combination (3 + 4>p eng /2. However, combining both measurements, it is possible in 
principle to extract all phases separately Thus (3 and 7 are determined and a can be solved 
for. Since all the measurements involve sine of some angle, there exists some ambiguity in 
determination of a definite angle. However the analysis involved in the process D°K* (892) 
is in principle expected to determine the definite value 7, if in addition one studies the 
exclusive processes Bd — > D°X° (X is K + tt~, K + tt~tt°, etc.) to remove discrete ambiguity 
I- 

We recall that in the SM with three generations, the sum of three CKM phases a, (3 and 7 
must be equal to tt. In order to check the validity of this unique feature, one would measure 
the CKM phases, for instance, through B& decay modes such as tttt, ipKs and D°K* (892) 
which are preferred experimentally and would yield a, (3 and 7, respectively, in the SM. 
However, as we can see from Table I, these modes would actually measure tt — (P+j + cj) Bd /2), 
P + ((pB d + 4>k)/2 and 7, respectively. The sum of these three angles would give tt + 4>k/2 
which can be a good indication for new physics unless 4>k turns out to be small. Even in 
case the experiments show the sum of these angles to be tt, there is still room left for extra 
physics because of the possible existence of 4>B d or 4>b s - Another interesting case is of multi- 
Higgs models, where SM phases might be absent. This corresponds to 7 = (3 = 0, a = tt. 
In that case, asymmetry in Bd —> ipK s is opposite in sign to Bd —>■ tttt, and 7 measurement 
will yield 0. 

If we concentrate just on the Bd decay modes, since these decay modes are more prefer- 
able from the experimental viewpoint, it is hard to extract all the CKM angles cleanly. But 
the angle 7 can still be measured without contamination of the extra phases. Since it seems 
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to be very difficult to extract a and (3 by using any independent methods, we suggest that 
a and (3 be determined using the unitarity triangle. Measuring the ratio of the CKM factors 
iKtdKi&l/lKdKfcl (e.g. by studying the spectra of charged leptons in the semileptonic pro- 
cesses b — > uu e e and b — > cu e e) and using 7 when measured, one can construct the unitarity 
triangle completely, which enables one to determine the phases a and (3 simultaneously. 
This angle (3 should be compared with the angle measured in the Bd decay modes such as 
ipKs in order to extract information about new physics. 

In conclusion, we have shown how the measurement of CKM phases as well as additional 
phases can be achieved when comparable contributions from beyond the standard model 
might be present. 

This work was supported by Department of Energy grant DE-FG03-96ER 40969. 
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TABLES 



TABLE I. B decay modes for measuring CP angles. 





Quark Process 


B d Modes 


Angles 


5 S Modes 


B s Angles 


(1) 


b — ► ccs 




(3+(4> Bd + <pK)/2 


DfD- 




(2) 


b — > sss 




P + {4>B d + <t>K + <Ppeng)/2 




(<pB s + <Ppeng)/2 


(3) 


b — > mZs 




— 


<M°, K+K- 




(4) 


6 — ► (ids 


7T K S , 




K°K° 


{<Pb s + <P' P en 9 )/ 2 


(5) 


6 — > ccd 


D+D-, Y>7r°, 
D°D° 


(3 + 0B d /2 


tpK s 


(0B S + <M/2 


(6) 


6 — > ssd 


K°K° 


{<f>B d + <en S )/2 


<j>K s 


+ (4>K + 4>B S + 4>peng)/2 


(7) 


6 — > uud, 
ddd 


7T7r, 7rp, 
7rai 


7r-(/3 + 7 + Bd /2) 


p°K s , tt°K s 


7 + {<t>B s + 4>k)/2 


(8) 


b — > cus, 
ucs 


D^P-ftT* (892) 


7 


D°cp<t> 
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FIG. 1. Plots of r(= AmB gluino / ArriB SM ) as a function of \i for different values of the gaugino 
masses (7711/2). The scalar mass mo=l TeV for the entire plot. 
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FIG. 2. Plots of r(= ArriB gluino / AniB SM ) as a function of the gaugino mass (mi/2) for different 
values of the scalar masses mo- 
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